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Abstract 

We examine the weakly interacting atoms in an ultracold Fermi gas leading to a state of macro- 
scopic coherence, from a theoretical perspective. It has been shown that this state can be described 
as a fermionic coherent state. These coherent states are the eigenstates of fermionic annihilation op- 
erators, the eigenvalues being anti-commuting numbers or Grassmann numbers. By exploiting the 
simple rules of Grassmann algebra and a close kinship between relations evaluated for more familiar 
bosonic fields and those for fermionic fields, we derive the thermodynamic limit, the spontaneous 
symmetry breaking and the quasi-particle spectrum of the fermionic system. 
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I. INTRODUCTION 



Dilute, ultracold, trapped atomic gases have served as model systems for exploring quan- 
tum phenomena at the fundamental level lH4|]. One of the basic issues is the condensation 
of a Fermi gas. Important experimental advances on alkali atoms 6 Li and 40 K J5, [f| have 
been made in this regard to study the interaction between the fermions. Since the inter- 
action between cold atoms characterized by s-wave scattering length can be controlled by 
magnetic field Feshbach resonances, it has been possible to realize both the repulsive (BEC) 
and the attractive (BCS) regimes and their crossover. In spite of a very large number of 



ytic solution of the many-body 
7H9|. However, the use of mean 



theoretical papers published over a decade, a complete ana 
problem along the BEC-BCS crossover still remains illusive 
field theory, pseudopotential and consideration of fluctuations around mean field has led 
interesting advancement in understanding the basic physics of ultracold atomic Fermi gases. 
For details we refer to the topical review by Giorgini, Pitaevskii and Stringari 

To capture the essential physics around the crossover it is necessary to comprehend the 
formation of Cooper pairs or molecule or BEC condensate or condensation of fermionic atom 

:he description of "Fermi conden- 



pairs 
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14j. A close look into this aspect reveals that t 



14j . This is because of the fact 



sate" around this crossover is still somewhat incomplete 
that a simple two-body physics of the resonance on the attractive side of interaction does 
not support a weakly bound molecular state - a point emphasized by Jin et. al. 12|. This 



implies that many-body effects must prevail in characterizing such a state. In the present 
study we focus on the description of this state in a weakly interacting Fermi gas. Pauli 
principle rules out the possibility of macroscopic occupation of a single quantum energy 
state of fermions. A question, however, remains, whether it is possible to realize a state 
of macro scop ic coherence among them. This is based on exploitation, following Cahill and 
Glauber [15] . of the close parallels between the families of the quasi-probability phase space 
density functions for bosonic and fermionic fields. A key element of the formulation is the 
fermionic coherent state defined as a displaced state where the displacement operator acts on 



the vacuum state. The transformation with the displacement operator disp 



aces a fermionic 



15Nl8] . Fermionic 



field operator over an anti-commuting number or Grassmann number 
coherent state is an eigenstate of annihilation operator, eigenvalues being the Grassmann 
numbers. This is reminiscent of the harmonic oscillator coherent state for which the dis- 
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placement operator displaces the bosonic field operator over a classical commuting number 



19 



21]. Our method for description of the macroscopic state of weakly interacting fermions 



is based on the fermionic coherent state and the associated algebra of the anti-commuting 



numbers, which have no clasical analog 15l-|l8|. They form the basis for demonstration of 
the thermodynamic limit, spontaneous symmetry breaking and the quasi-particle spectrum 
of the atoms comprising the coherent state. A major goal of the approach is to show that it 
is possible to achieve a coherent state of weakly interacting fermions for which the scatter- 
ing length can be either positive or negative. And the state can exist as independent entity 
irrespective of the details of Cooper pairing or molecule formation. 

The paper is organized as follows: In Sec. II we briefly introduce the basic model for 
weakly interacting many body systems under ultracold condition. Since Grassmann num- 
bers play a crucial role in the formulation of the problem we briefly review in Sec. Ill the 
relevant aspects of Grassmann algebra before introducing the fermionic coherent state fol- 
lowing Cahill and Glauber flo| . That the state of macroscopic coherence can be realized 
as a fermionic coherent state is the main theme of this section. We further examine the 
aspect of thermodynamic limit and spontaneous symmetry breaking . In Sec. IV we derive 
the quasi-particle spectrum using Bogoliubov transformation. The paper is concluded in 
Sec.V. 



II. WEAKLY INTERACTING FERMI GAS 

To begin with we consider a dilute gas of fermionic atoms. The interaction between them 
is as usual nonzero only when the range of inter-particle interaction ro is much smaller than 
the average inter-particle distance d, i.e., r << d = p -1 / 3 = (y) 1 ^ 3 - This also ensures that 
the properties of the system can be expressed in terms of a single parameter, scattering 
length, a. The diluteness criterion in terms s-wave scattering length can be expressed as 
la)/? 1 / 3 << 1. We also assume that the temperature of the dilute gas is so low that the 
momentum q as distributed thermally is much smaller than the characteristic momentum 
q c = h/r . The scattering amplitude is independent of momentum q when q << q c . While 
the temperature for onset of quantum degeneracy goes as p 2//3 , temperature for BEC-BCS 
crossover is of the order of Fermi temperature for a range of dimensionless fcp|o| values, 
where k-p corresponds to Fermi wave vector j?-9|. 
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The Hamiltonian of the system can be expressed in terms of the field operator ip(r) as 
H= / (— V$ (r)Vip(r)) dr + - I ^\r)$ (r')V {r - r')^{r)${r')dr dr' (2.1) 



\2m * v ' j 2 

where V(r — r') is the two-body potential. For the moment we have not included the 
external fields. In absence of external fields, the iV particles move only as a result of their 
mutual interactions. The field operator V>(r) annihilates a particle at a position r and can 
be expressed as ip(r) = ^ 0j(r)aj. ci-j(aj) is the annihilation (creation) operator of a particle 
in the single particle state <pi{r) and these operators obey anti-commutation relations 

{Si, a]} = 8ij] {Si, cij} = {a\, a]} = (2.2) 

The wave functions 0j(r) satisfy orthonormal condition 

#(r)^(r) = 5 is (2.3) 

The field operator then follows anti-commutation relation {ip(r), fo(r')} = ^ 0i(r)0* (r') 
= <5(r — r'). Now expressing the field operator V>(r) and the interaction term, V(r,r') (= 
V(r — r')) as ip(r) = J2 k a-k^^r and V g = J V(r)e~~*~dr, we may simplify the Hamiltonian 

as 

h 2 k 2 1 

H = X) "^T^ + ^7 Yl V A+A2- q ^ 2 a kl (2.4) 

k fci ,&2 ,<7 

We take into account of the momentum conservation during a two-fermion interaction 
such that every single term within the summation describes the annihilation of a pair of 
fermionic atoms and the creation of another pair. Since the description of the macroscopic 
properties of the gas does not require the detailed form of the two-body interaction, mi- 
croscopic potential V(r) may be replaced by a soft, effective potential V e ff(q) which, in 
turn, gives the low momentum value of its Fourier transform Vo- For such small momenta 
q « q c = h/r we are allowed to consider Vq or the s-wave scattering length 'a' becomes 
independent of its momentum q |z-2|. 

The Hamiltonian therefore reduces to 

H = 1^ ~2^ a * ak + ^7 L < 1+q < 2 - q ak 2 a kl (2.5) 

k ki,k 2 ,q 

This is the starting Hamiltonian for weakly interacting Fermi atoms and henceforth will be 
used in the following sections. Finally we emphasize that the Hamiltonian considered here 



( Eq , 2.1-2.5) follows Greiner |21| and does no. refer to speeine spin states. Th.s is an 
important point of departure from many theoretical approaches to quantum Fermi liquid [7fl 
and BCS-BEC crossover jj-fl. 

To proceed further for a systematic description we first separate out the field operator 
V>(r) into the coherent term (i = F) and the incoherent term (i ^ F) as; 



i(r) = <t) F (r)a F + ^ <M 



r a; 



(2.6) 



22] for bosons and 



We do this separation by taking a hint from Bogoliubov approximation 
the close parallelism between the expressions developed for bosonic and fermionic fields by 
Cahill and Glauber [15] . We may note that an implementation of Bogoliubov approximation 
J22I for fermionic fields requires replacement of ap and a F by anti-commuting numbers, say 



y-p and y^ |15Hl8|. One can always multiply these numbers by a numerical phase factor 
as y-pc 16 and y^,e~ %9 without changing any physical property This phase plays, as we 
will see, a major role in characterizing the coherence in the weakly interacting ultracold 
Fermi gas. This reflects the gauge symmetry exhibited by all the physical equations of the 
problem. Making an explicit choice for the value of the phase actually corresponds to a 
formal breaking of gauge symmetry. 

Equivalently the Bogoliubov approximation for fermions is equivalent to treating the 
macroscopic component 0F(r)aF as a "classical" non-commuting field ij}-p{r) so that ip(r) 
may be rewritten as 



(2.7) 



The above ansatz for the fermionic field operator can be interpreted as the expectation 
value (V>(r)) different from zero. This is not possible if the coherent state is a particle 
number eigenstate. Although such coherent states are well-known for bosonic fields and 
have been the basis for understanding BEC 23N27|. an extension of the scheme to their 



fermionic counterpart is not straight-forward. The primary reason may be traced to a basic 
issue. Since fermions anti-commute their eigenvalues must be anti-commuting numbers, 
as pointed out by Schwinger [16]. Such numbers are Grassmann numbers which can be 



dealt with by simple techniques of Grassmann algebra 



171 . Il8| . In what follows in the next 



section we make a little digression on this algebra centering around fermionic coherent state 



following Cahill and Glauber 



15] which forms the integral part of the description of the 



interacting fermions near the crossover. 

Although an important guideline of the present formulation is the close parallels between 
some aspects of bosons and fermions it is important to highlight the other essentials. The 
following digression may be interesting. While for a bosonic gas the condensed phase is 
formed by the particles in the zero energy ground state, in the fermionic gas the particles 



may "condense" [28 
Kothari and Nath 



in the Fermi level. Such a prediction was made many decades ago by 
28j in course of examination of Born's reciprocity principle [29]. In the 



present scenario, we note that the BEC-BCS crossover temperature for a Fermi gas is less 



than or equal to the Fermi temperature 



1CH14|. a point that hints towards this assertion. In 



what follows we will show that a macroscopic state of weakly interacting ultracold fermions 
can give rise to a fermionic coherent state. 

III. GRASSMANN ALGEBRA AND FERMIONIC COHERENT STATES; A CON- 
NECTION TO MACROSCOPIC COHERENCE OF FERMIONS: 

A. Grassmann variables and their properties: 

We now summarize some of the properties of the anti-commuting classical variables rele- 



vant for our future discussions. These variables can be treated wit 



lin the scope of Grassmann 
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They possess very un- 



algebra, which are well studied in mathematics and field theory 
common properties. 

Let, y = {yi},i = 1,2, n. define a set of generators which satisfy anti-commutation 

properties, 

ViUj + VjVi = {V%, Vj} = (3- 1 ) 

This, in particular, implies that yf = 0, for any given H\ Since the square of every Grass- 
mann monomial vanishes, a nonzero Grassmann monomial can not be an ordinary real, 
imaginary, or complex number. In other words, Grassmann variables are nilpotent, an im- 
portant property for the treatment of fermions The anti-commuting numbers yi and 
their complex conjugates y* are independent numbers and satisfy 

buy*} = (3.2) 
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They also anti-commute with their fermionic operators 

{y i} a j } = 0; {y h a]} = (3.3) 

And hermitian conjugation reverses the order of all fermionic quantities, both the operators 
and the anti-commuting numbers. For instance, we have 

(amalvl^ = V4a 3 y*a\ (3.4) 

An analytic function of only one Grassmann variable can be expressed as a simple Taylor 
expansion 

f(y) = fo + fiV (since y 2 = 0) (3.5) 
Thus exponential function for Grassmann variables has only two terms 

OO 

eV = J2li =1+ y ( since ?/ 2 = o) (3-6) 

fc=0 



Therefore for a function of single Grassmann variable, integration (integration due to Berezin 



17j) is identical to differentiation. The fundamental rules of integration over the complex 



Grassmann variables are as follows; 

r dy l = ldy* = (3.7) 

Sa (3.8) 



J dyi yj = 6ij J dy* y* 



This difference between ordinary variables and Grassmann variables has many interesting 
consequences. We are typically concerned with pairs of anti-commuting variables yi and y*, 
and for such pairs we will confine ourselves to the notation 



d yi = j dy*dyi (3.9) 
in which the differential of the conjugated variable dy* comes first and we keep in mind that 

dy-dyi = -dyidy* (3.10) 
We can write the multiple integrals over such sets as 

j d 2 y = j Y[dy*d yi (3.11) 
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B. Fermionic coherent states: 



We are now in a position to introduce the fermionic coherent states. In analogy to har- 
monic oscillator coherent state \ct) 19| defined as a displaced state where the displacement 
operator D(a) = exp ^X/j( a i^l ~ a j^i|j ac t s on the vacuum |0) as \oc) = D(a)|0), 
being a set of complex numbers 19M2l|. it is possible to construct a displacement operator 
for fermions |15| as 



D(y) = exp (^(alvi - y*di) j 



(3.12) 



for a set of y = {y{\ Grassmann variables and the fermionic coherent state can then be 
constructed by the action of this displacement operator (Eq. 3.12) on the vacuum state as 

|y) =D(y)|0) (3.13) 



1. Displacement operator and its properties: 



An important property of the Grassmann variables is that on multiplication with 



fermionic annihilation or creation operators, their anti-commutivity cance 
erators. Thus the operators a\yi and y*dj simply commute for i ^ j 
displacement operator may be written as the product 



s that of the op- 



15]. Therefore the 



D(y) = JJexp (a\yi - y*a^j 

i 

= n 



1 + a\yi - y*a,i + \a\ai - - j y*y l 



(3.14) 



Similarly the annihilation operators at commutes with all the operators a\yi and y*di for 
k 7^ i, and so we may compute the displaced annihilation operator by ignoring all modes 
except the k-th one 

D t (y)«fcD(y) = JJexp (y*di - a\y^j .a k . JJexp (afy - y*hj 



exp i y* k a k - d[y k ) .d k . exp (d[y k - y* k a k 



ak + yk 



(3.15) 



Similarly we have 



6^)4%) = a\ + yl 



(3.16) 
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2. Properties of the coherent states: 



By using the displacement relation [Eq (3.15)] we may show that the coherent state is an 
eigenstate of every annihilation operator a^; 

a fc |y) = a fe D(y)|0) = D(y)D t (y)a fc D(y)|0) 
= D(y)(a jfc + y fe )|0) 
= y fc D(y)|0) = 2/*|y> 



(3.17) 



The adjoint of the coherent state |y) can be similarly defined as (y\a k = (y\yl- 
The inner product of two coherent states is 



(y'|y> = exp (v^yfyi - + vtViYj 



(3.18) 



and using the completeness properties of the coherent states, any arbitrary coherent state 
|y') can be expanded 15|, [l8j as 



|y') = J rf 2 y(y|y'>|y> 



(3.19) 



which immediately follows from the resolution of identity 



J rf 2 y |y)(y|=i 



(3.20) 



C. Weakly interacting Fermi gas in coherent state representation: 

The key idea behind the present formulation is whether a state of macroscopic coherence 
of an weakly interacting Fermi gas can be achieved as a fermionic coherent state. Such 



an approach gives the result which is equivalent to that of t 
of particles iV in the limit N — > oo for the BEC case 123 



re problem of fixed number 



27fl. In the same spirit, we 



■> oo for the BEC case 
can extend the coherent state approach to its fermionic counterpart. Before proceeding 
further, we discuss the physical states and operators, which will be helpful for the future 
development. 
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1. Physical states and operators 



According to Cahill and Glauber, 15] a state is physical if it changes at most by a 



phase when subjected to a rotation of angle 2n about any axis, 

U(n,27r)|-0) = e ie \^) (3.21) 

Since, fermions carry half-integer spin, a state of any odd number of fermions changes 
by the phase factor -1, while even number of fermions are invariant under such 27T rotation. 
Thus physical states are linear combinations of states with even number of fermions or linear 
combination of odd number of fermions. On the other hand, a linear combination of states 
|0) and |1) like ^(|0) + |1)) must be unphysical since it changes to a different state under 
27r rotation. 

An operator is physical if it maps physical states onto physical states. Physical operators 
are either even or odd. In all physical contexts, if we consider = Ylk ^l^fc as the total 
fermion number which is conserved, then any state arising from the eigenstate of N must 
remain an eigenstate of N. This can be derived on the basis of U(l) invariance of all the 
interactions under the transformation U(9) = e t9N . We have the transformation on fij and 
a\ as follows; 

e -MN he i9N = ^ (3>22) 

and 

e -i6N & yN = e -ie a t (3.23) 

for fermion conserving interactions that involve and a\; and the phase factors get can- 
cel out. We now emphasize that the coherent states undergo a simple change under this 
transformation 

^)|y) = e^|y) = |e 4 V), (3.24) 
However the scalar product remains invariant 

(e 4 VKV> = <y|y> (3.25) 
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2. Spontaneous symmetry breaking: 

We now note that one can always multiply the coherent state by an arbitrary phase 
factor e 10 without changing any physical property (as far as physical states are concerned). 
This is the manifestation of gauge symmetry in the problem. Physically, the lack of a force 
responsible for phase stabilization of the system is the origin for the random phase of a 
condensate. However, when one refers to Fermi systems the macroscopic system is expected 
to choose spontaneously a particular phase Making an explicit choice for the phase 
in spite of the lack of a preferred phase value (referred to as a spontaneous breaking of gauge 
symmetry) implies that the macroscopic state is in or close to a coherent state. From the 
symmetry point of view, the situation is quite interesting and can be further elaborated as 
follows. 

The states |y) are not invariant under the number operator N = X^fc^l^fe) w hile the 
Hamiltonian (Eq. 2.5) commutes with N, i.e., 

e *#|y) = | e 4 V); e M He" M = H (3.26) 

The operator e t9N applied to |y) produces a state with the same energy, with a phase shifted 
by u 9". Since the overlap of coherent states obey (y'|y) = exp (J2iV'i*yi ~ \{y'iV'i +?/*?/*)) 
[Eq. (3.18)], any two different states |y), |y') with different phase factors, |y), |e j6, y), param- 
eterized by a phase variable < 9 < 2tt, are macroscopically distinct. While the microscopic 
Hamiltonian (H) [Eq. 2.5] has global U(l) symmetry, the state does not possess such sym- 
metry, since adding a phase factor to the state |y) produces a different state altogether. 



3. Population fluctuation and phase stabilization: 

One can easily calculate the population fluctuation in the coherent state given by Eq. 
(3.13) as 

(AN 2 ) = (N 2 ) - (iV> 2 = £ VkVh = (N) (3.27) 

k 

Here we have to use the property [15] of the fermionic number operator N = a\ak and 
the nilpotency of the Grassmann variables i.e., (y^Uk) 2 = 0, for each k. 
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Then the phase fluctuation in the coherent state is given by 

(Ad 2 ) = K- = -L- = -— ^ (3.28) 

4(AiV2) 4(iV) tEkVkVk 

Since N and 8 are the canonically conjugate pairs, the phase is stabilized (AO 2 ) « 1, only 
at the cost of the enhanced fluctuation in population i.e., (AN 2 ) » 1. Since the coherent 
state can be expanded as a coherent superposition of the particle number eigenstates [Eq. 
(3.13)], the constructive and destructive interferences among different number eigenstates 
result in the stabilized phase but with finite particle number noise. Population fluctuation 
and phase stabilization are the typical signatures of spontaneous symmetry breaking 27]. 



4- Thermodynamic limit; Grassmann-Bogoliubov approximation: 

Once a macroscopic state of interacting fermionic atoms is realized as the coherent state, 
neither a F nor d F annihilates the state. Since, we are interested in the behaviour of a gas 
of fermionic atoms, i.e., in the large particle number and volume, it is necessary to consider 
the so-called thermodynamic limit N — > oo, V — > oo but with constant density p = y. 
In this limit, the anti-commutation relation between fermionic operators a F and a F becomes 

{of^} = a F a F + a F a F = \_ _^ Q ^ y _^ ^ ^ 

In the limit, V — >■ oo, we are allowed to forget the operator character of ap and a F 
and they can be replaced by numbers, i.e., we obtain the "classical" limit of the fermionic 
operators. To make this point explicit, let us define 

a F a F = N F (3.30) 

where N F represents the number operator for the particles in the coherent state. Because of 
(y|^F|y) = UfUf 7^ and the anti-commuting properties of Grassmann variables (Eq. 3.1), 
it follows 

(y|{a F ,4}|y) _ 1 + {2/f,S/f> 1 ^ n /o qi \ 

v " v — v {6 - 6L) 

This leads to the natural starting point of what we may call the Grassmann-Bogoliubov 
approximation. Here one can identify y ¥ y-p = N-p(m N) as the average particle number of 
the state of macroscopic coherence in the thermodynamic limit. 
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One point is to be noted here. One may conclude that ant i- commutation obeyed by 
fermions do not have a classical analogy since they do not fulfill classical Poisson brackets 
as obeyed by bosonic commutation relations. But this should not lead to misunderstand- 
ing. The number operator N = Nk and the Hamiltonian operator H = Y2k e k^k have 
classical limits because they are bilinear in d&, at. For the particles obeying Fermi-Dirac 
statistics, only quantities like charge, energy, current density or number density are measur- 



able classically because they are bilinear combination of field amplitudes d^, a k 18|. The 
amplitudes of the Fermi field is linear in dk, dj. and can not be measured classically. Ex- 
trapolating the idea a bit further we can emphasize that Grassmann fields themselves and 
fermionic field operators are, by construction, fermionic while a product of even number 
of fermionic quantities or Grassmann variables is bosonic which makes it experimentally 



relevant 
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18 



30|. 



With this, we now return to our starting Hamiltonian (Eq. 2.5) characterizing the weakly 
interacting Fermi atoms. In the next section the energy spectrum is calculated order by order 
under Grassmann-Bogoliubov transformation. 

IV. QUASI-PARTICLE SPECTRUM AND QUANTUM FLUCTUATIONS: 
A. Lowest-order approximation: 

In the first approximation, we can neglect all the terms in the Hamiltonian (Eq. 2.5) 
containing the operators d^ and a\ with k ^ F. This implies that the quantum fluctuation 
(5ip(r)) in Eq. (2.5) can be ignored. Under such approximation the coherent state must 
be occupied by macroscopically large number of particles i.e. N-p ~ N. The replacement 
of dp and d F by '?/fC 1 ^' and '2/ F e~ ieh respectively, then becomes quite straight forward. 
This substitution can not be made for a realistic potential since it would result in a poor 
approximation at short distances of order r , where the potential is strong and correlations 
are important. The replacement is instead accurate in the case of a soft potential whose 



perturbation is small at all distances 
order therefore takes the form: 



27j . The energy of the interacting system in the lowest 



= ~2m ? a¥ + ^7 a F a F a Fa F = ~2^y-pV^ + (4.1) 
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To the same order of approximation, one can easily express the parameter V of Eq (4.1) 
in terms of the scattering length a, using the result Vo = 4n ^ a , under Born approximation. 
However, due to the anti-commuting nature of the Grassmann variables [Eq. (3.1)] the 

fc2 £.2 

energy is given by only the first term in Eq. (4.1), i.e., E F = -^-^N-p — ^N F , where, 
= VyVf represents the average number of particles and is the chemical potential of 
the weakly interacting fremions in the coherent state. 



B. Higher order approximation; calculation of quasi-particle spectrum: 

The result E F = /j,N F for the energy of the interacting fermions in the lowest order has 
been obtained by taking into account of Eq (2.5) only for the particle operators a k and a\ 
with k — F. The terms containing only one particle operator with fc ^ F do not enter 
into the Hamiltonian (Eq. 2.5) because of momentum conservation. By retaining all the 
quadratic terms in the particle operators with k ^ F, in the next higher order, we obtain 
the following decomposition of the Hamiltonian; 

T'T \ V At A ^0 A f At -S VO V— V / t A t„ , a a "t "t \ , S 

' k(Xk + 2^*^ + 2V 2-*t (^ 4a F a fc a F«fc + a ] k a ] _ k a F a F + a F a F a k a^ k J (4.2) 

k k^F 

Now, the replacement of a F and a F by y F e~ %e and y F e ld , as carried out previously, yields 
the following expression for the Hamiltonian; 

H = Yl Y^^ k + ^ y v y v yFyF + ^7 Yl ( 4 ^F«IyFa fe + a[a<_ k y F y F + y F y F a k a^(A.3) 

k fc^F 

The third term of this Hamiltonian represents the self-energy of the excited states due 
to the interaction, simultaneous creation of the excited states at momenta k and —k and 
simultaneous annihilation of the excited states, respectively. But, the simultaneous creation 
and annihilation of the excited states do not contribute to the Hamiltonian (H) due to the 
Grassmannian nature of the fermionic field variables. Introducing the relevant interaction 
coupling constant g fixed by the s-wave scattering length a as 

4ich 2 a 

9 = 4.4 

m 

the Hamiltonian (Eq. 4.3) can be written as 

H = E F + \j^(E k + 2pg)(a[a k + al k a. k ) (4.5) 
k 
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Here s *S n indicates that the terms k = F are omitted from the summation, p represents 
the fermionic density expressed as p = = y). The physical coupling constant g 

renormalizes the effective potential Vo- E k refers to the energy when the interaction g = 
and is given by Ef. = Unlike the bosonic case, the Hamiltonian in Eq (4.5) is peculiar 
in a sense that it has no terms with creation operators only or annihilation operators only. 
In other words, bosonic Hamiltonian does not conserve the number of particles. On the 
other hand in all physical contexts that have been explored experimentally, the number of 
fermions or more generally, the number of fermions minus the number of antifermions is 



strictly conserved 15]. This conservation law leads to further restriction on the permissible 
states of the field. If a system starts with a state of fixed number of fermions, the conserva- 
tion law restricts the set of accessible states considerably more than the 2tc superselection 
rule mentioned earlier. Transitions can not be made, for example with different even number 
of fermions or between states with different odd number of fermions. 

We now look for a solution of the problem (Eq. 4.5), i.e., energy eigenvalues of the Hamil- 
tonian. Since the coherent state, represents a combination of unperturbed eigenfunctions, 
neither a k nor a\ annihilate this state. The problem can be solved exactly by a canonical 
transformation, namely Grassmann-Bogoliubov transformation. To this end we introduce 
new operators 

a k = u k A k + v*_ k A ] _ k 

4 = u l A l + v ~kA. k (4.6) 

This transformations introduces a new set of operators A k and A k on which we impose the 
same fermionic anti-commutation relation 

{A k ,Al} = 6 ktk , (4.7) 

as obeyed by the original particle operators a k and a k , to make the transformation canonical. 
It can be easily verified that Eq. (4.7) are satisfied if 

\u k \ 2 + |f_ fc | 2 = 1 (for any k) (4.8) 

Therefore u k and v_ k can be chosen parametrically as 

u k = cos8 k and v_ k = sin9 k (4.9) 
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By inserting Eq (4.6) into Eq. (4.5) we obtain 

H = E F + \Y!( Ek + 2 P9) 2 \ v -k\ 2 + \Yl( Ek + 2 P9)-(\ u k\ 2 - \v-k\ 2 )-{A\A k + Al k A- k ) 
k k 

+\Y!( Ek + 2 P g)2u k v„ k {A\Al k + A k A_ k ) (4.10) 
fc 

In order to make H diagonal in A k and A\ we use the freedom to eliminate the last term of 
Eq. (4.10) i.e., 

\Y^^ Ek + 2 P9) 2 u k V-k = (4.11) 

k 

With Eq. (4.9) and by defining \(E k + 2pg) as ctfe, Eq (4.11) gives for each k-th mode 

a fe sin20 fe = O (4.12) 

The above condition is satisfied both for positive to negative values of a k . The allowed 
values of the coefficient 9 k are 

TfiTl 

6 k = ±— (m=integer) (4.13) 

It is easy to note that for the above 9 k values, the first term of Eq (4.10) corresponding to 
E we obtain 

E = E F + ^( E k + 2 P9)^k\ 2 

k 

= E F + \j^(E k + 2pg)2 sin 2 26 k (4.14) 

k 

Since \Y^! k {E k + 2pg)2 sin 2 29 k = we have E = E F . This is the same result in the lowest 
order approximation obtained earlier. So the H in Eq. (4.10) is modified as 

H = E F + \j^(E k + 2pg).(\u k \ 2 - \v_ k \ 2 ).{A\A k + Al k A_ k ) (4.15) 

k 

Identifying \{E k + 2pg).{\u k \ 2 — \v^ k \ 2 ) = a k cos2# fe = a k , the Hamiltonian in its final form 
can be written as 

H = E F + lj^(E k + 2pg)(AlA k + Al k A_ k ) 



2- 

k 



= E F 4- ^e k (A[A k + Al k A_ k ) (4.16) 



2- 

k 



16 



where the quasi-particle energy e k is given by 

e k = ^ + 2pg = E k + 2pg (4.17) 
2m 

E k = 4^ stands for the free particle energy if the interaction g = 0. The operator A^A/- 
resembles a particle number operator and has eigenvalues and 1. Hence the coherent state 
is determined by the requirement that 

A fc |y(0)) = O for all A; ^F (4.18) 

Furthermore, all quasi-particle states correspond to different numbers of non-interacting 
fermions. The expression (4.16) is the fermionic counterpart of the Bogoliubov Hamiltonian 
for bosons. E-p term originating purely as kinetic energy term for bosonic case is zero since 
the coherent state corresponds to k = (instead of k = F). Secondly the ground state is the 
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vacuum of the quasiparticle operators A k , A k . In general, the spectrum is gapped 
everywhere corresponding to a non-zero difference between the first excited state and the 
ground state. The theory does not put any restriction upon the sign of g and hence it can 



be both positive and negative 
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U-ll4j. For a critical negative g i.e., g c it is possible that e k 
becomes zero when E k = —2pg c . In this case the spectrum can be gapless [30]. This is in 
consistent with the general perception that in the case of spontaneously broken continuous 
symmetry, there is always a low-lying excitation energy e k which satisfies e k — > as its 
momentum k — > 0. In the present context quasi-particle energy e k smoothly goes to 
zero without a gap as the momentum goes to zero, even though the dispersion relation is 
quadratic in momentum rather than linear unlike the Bose gas. 

The coherent state \y(9)) is annihilated by all A k , (k ^ F). The transformation (4.6) 
can be represented as 

A k = #(#K#(#)t (4.19) 
A\ = ti(6)alti(6y (4.20) 

[&) represents the unitary displacement operator which produces the coherent state as 
given in Eq (3.13). 

|y(0)) = tf(0)|O) (4.21) 
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The connection between tyf(6) and D(y) in Eq. (3.13) can be established as follows. 

#(0)|O) = #(0).D(y)|O) 

= u(0)\y) 

= exp (iONj |y) 

= \e ie y) (4.22) 

Under the rotation of all Grassmann variables y^ by the same angle i.e., the phase 
transformation of y^ — > e %6 yk 



15] the displacement operator produces altogether a different 
coherent state \e %e y). But for a physical rotation is allowed to run from to 2tt. So, the 
permissible values of '6^ for which Eqs (4.19) and (4.20) are valid are 9 = 0, ±~, ±7r, ±^ as 
obtained from Eq. (4.13). These values of l 9' are the same as obtained earlier by Braungardt 
et. al. in the context of Fermi gas in optical lattice in Id ^30 1 . For these specific choices of 
the phase of the coherent state gets stabilized. 

V. CONCLUSION 

In this article we have tried to understand the coherent state of a weakly interacting 
ultracold Fermi er close analogue of Bose condensate. Notwithstanding statistical 

differences, the coherent state of fermions bear a close kinship with that for bosons. Since 
the fermionic operators anticommute their eigenvalues are anticommuting numbers. These 
numbers or Grassmann variables play an important part in the formulation of the present 
theoretical scheme. The fermionic coherent state is not invariant under rotation e l9N with 
number operator N for fermions, while the Hamiltonian remains unitarily invariant. For pre- 
ferred choices of 8, the coherent and the rotated coherent state are therefore macro scopically 
distinct states because of spontaneous symmetry breaking. The description of the coherent 
state of weakly interacting fermions and the rules of Grassmann algebra allow us to consider 
the appropriate thermodynamic limit of the system. Bogoliubov approximation can be im- 
plemented within the framework of this scheme to realize this coherent state as a state of 
macroscopic coherence of fermions. Unlike the bosonic case the lowest order approximation 
on the Hamiltonian describing the weakly interacting fermions yields an energy contribution 
due to kinetic energy of the particles only. The quasi-particle spectrum arising from the 
higher order interaction exhibits in general, the gaps Our analysis reveals that 
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the close parallels between the phase space quasi-probabilities of the boson and fermions 
as pointed out by Cahill and Glauber [15| have deep rooted consequences in the physics of 
ultracold degenarate Fermi gas. 
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